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We perform an indirect search for Weakly Interacting Massive Particles (WIMPs) using the 
MACRO detector to look for neutrino-induced upward-going muons resulting from the annihilation 
of WIMPs trapped in the Sun and Earth. The search is conducted in various angular cones centered 
on the Sun and Earth to accommodate a range of WIMP masses. No significant excess over the 
background from atmospheric neutrinos is seen and limits are placed on the upward-going muon 
fluxes from Sun and Earth. These limits are used to constrain neutralino particle parameters from 
supersymmetric theory, including those suggested by recent results from DAMA/Nal. 

95.35, 14.80.-j 



I. INTRODUCTION 

There are many hints for the existence of non-baryonic dark matter in our universe [jjj , which may consist of Weakly 
Interacting Massive Particles (WIMPs). Measurements of the mass density of the universe, Qm, indicate a density 
considerably in excess of the density of baryonic matter allowed by Big Bang nucleosynthesis (BBN) calculations. 
For example, the POTENT analysis of peculiar velocities of galaxies Q excludes CIm < 0.3 at the 2.4<r level (Om is 
in units of the critical density). In contrast, the limits for baryonic mass (fib) from BBN are 0.005 < < 0.10 at 
95% c.l. H for 0.4 < h < 1, where h is the scaled value of the Hubble constant. This discrepancy is an indication 
that there must be some undiscovered non-baryonic dark matter. Large-scale structure models coupled with COBE 
data Q also favor cosmological scenarios with large amounts of dark matter. One recent model favors a mixture of 
70% Cold Dark Matter, 20% Hot Dark Matter and 10% baryons J|. At smaller scales (roughly 100 to 10000 kpc) 
virial estimates on groups and clusters of galaxies and rotation curves of spiral galaxies || are consistent with smaller- 
values of SIm- Furthermore, dark matter in galaxies is also motivated by the fact that dark matter halos seem to help 
stabilize spiral disk structure Q. Nevertheless, the amount of matter indicated cannot be accounted for by baryonic 
matter: MACHOs are unable to account for all of the dark matter halo of our galaxy (7). 

A long list of Cold Dark Matter nonbaryonic candidates has been suggested, among which the supersymmetric 
(SUSY) neutralino, considered in this paper, and the axion seem to be the most promising ||. SUSY postulates a 
symmetry between bosons and fermions predicting SUSY partners to all known particles (for a review see Ref. Q). 
SUSY solves a host of particle physics questions such as the "hierarchy problem" (explaining the large difference 
between the weak and GUT scales); generating electroweak symmetry breaking through the Higgs mechanism; and 
stabilizing the Higgs mass at the weak scale. In theories where R-parity is conserved there exists a stable Lightest 
Supersymmetric Particle (LSP). If the neutralino is the LSP it is a natural WIMP candidate: it is a weakly interacting 
particle with a mass between roughly a GeV and a TeV and would be expected to have a significant relic density. 

The neutralino x is the lightest linear superposition of gaugino and higgsino eigenstates (gauginos and higgsinos 
are SUSY counterparts to the Higgs and gauge bosons): 

X = 017 + a 2 Z + a 3 Hi + a^H 2 (1-1) 

where 7 and Z are gaugino states; H\ and H2 are higgsino states. If the minimal supersymmetric extension of the 
standard model (MSSM) ^ is coupled with some GUT assumptions (Mi = 5/3M 2 sin 2 9w where 0w is the Weinberg 
angle), then the neutralino mass depends on one of the two gaugino mass parameters, M\ or M2, on the higgsino 
mass parameter /1, and on the ratio of the Higgs doublet vacuum expectation values, tan/3. The phenomenology of 
neutralinos is determined by the composition parameter P — a\ + a\. Some other parameters must be determined 
in order to define the processes induced by neutralinos, because Higgs and the supersymmetric partners of fermions 
(sfermions) play a relevant role. In the MSSM there are two Higgs doublets, hence three neutral Higgs fields (two 
scalar and one pseudoscalar) . The Higgs sector is determined by two independent parameters: tan/3 and tcia, the 
mass of the pseudoscalar neutral boson. The other parameters in the Lagrangian of the model are the bilinear and 
trilinear parameters connected to the spontaneous symmetry breaking. The number of parameters needed to describe 
neutralino phenomenology may be reduced further by assuming that all the trilinear parameters are zero except for 
the third family, which are assumed to have the common value, A, and that all the squarks and sleptons are degenerate 
with common mass, mo- In this paper we use a model developed following the assumptions above p0[ . 

Currently, supersymmetric parameter space is constrained by accelerator searches. Data from LEP |llj lead to 
lower limits on the neutralino mass, m x , between 20-30 GeV. These limits are model-dependent and correlated to 
limits from chargino (a mixture of w-inos and charged higgsinos) searches. 

In this framework, "direct" and "indirect" methods for detecting Galactic halo WIMPs can probe complementary 
regions of the supersymmetric parameter space, even when more extensive LEP 2 results become available. Direct 
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methods detect WIMPs via a direct interaction of a WIMP, such as by observing the energy deposited in a low- 
background detector (e.g., semiconductors or scintillators) when a WIMP elastically scatters from a nucleus. Indirect 
methods look for by-products of WIMP decay or annihilation such as neutrinos resulting from the annihilation 
of WIMPs. An excellent prospect for indirect WIMP searches is to look for high energy neutrinos from WIMP 
annihilation in the core of the Earth or the Sun In MACRO, such neutr inos would be detected as neutrino- 

induced upward-going muons which can be distinguished from downward-going cosmic-ray shower muons. 

II. INDIRECT SEARCH THROUGH UPWARD-GOING MUONS 

Dark matter WIMPs in the Galactic halo can be captured in a celestial body by losing energy through elastic 
collisions and becoming gravitationally trapped. As the WIMP density increases in the core of the body, the WIMP 
annihilation rate increases until equilibrium is achieved between capture and annihilation. High energy neutrinos are 
produced via the hadronization and decay of the annihilation products (mostly fermion-antifermion pairs, weak and 
Higgs bosons) and may be detected as upward-going muons in underground detectors. 

There are many calculations of the expected neutrino fluxes from WIMP capture and annihilation in the Sun and 
Earth || |l5|Jl6f| . The capture rate for an astrophysical body depends on several factors: the WIMP mean halo velocity 
(~ 270 km s _1 ); the WIMP local density (p x ~ 0.3 - 0.6 GeV cm -3 ); the WIMP scattering cross-section; and the 
mass and escape velocity of the celestial body ||. The WIMP may scatter from nuclei with spin (e.g. hydrogen in 
the Sun) via an axial- vector ( "spin-dependent" ) interaction in which the WIMP couples to the spin of the nucleus or 
via a scalar interaction in which the WIMP couples to the nuclear mass. In axial-vector interactions the probability 
for a given energy loss is constant up to the kinematic limit of the interaction, while for a scalar interaction there will 
be a suppression of the cross section at high momentum transfers ||. Elastic scattering is most efficient when the 
mass of the WIMP is similar to the mass of the scattered nucleus. Hence, the heavy nuclei in the Earth make it very 
efficient in capturing WIMPs with m x < 100 GeV (the resonance effect |l^]). The Sun, in contrast, has a smaller 
average nuclear mass, but is nonetheless efficient in capturing WIMPs due to its larger escape velocity. 

A WIMP annihilation signal would appear as a statistically significant excess of upward-going muon events from 
the direction of the Sun or of the Earth among the background of atmospheric neutrino- induced upward-going muons. 
High energy neutrino-induced upward-going muons tend to retain the directionality of the parent neutrino. This 
directionality permits a restriction of the search for WIMP annihilation neutrinos to a narrow cone pointing from 
the Earth or Sun, greatly reducing the background from atmospheric neutrinos. This effect, along with the increase 
in neutrino cross-section with energy and longer range of high energy muons, means that this method of detection 
achieves an increasingly better signal to noise ratio for high WIMP masses. 

Data on upward-going muons from the core of the Earth and of the Sun have been presented by several experiments, 
notably Baksan JTsfl , Kamiokande [[t9| and 1MB pOfl . In this paper we describe a WIMP search using the MACRO 
detector. 



III. WIMP SEARCH IN MACRO FROM THE EARTH AND THE SUN 

The MACRO apparatus [|T) , located in the Gran Sasso Underground Laboratory of the Italian Istituto Nazionale di 
Fisica Nucleare, detects upward-going muons using a system of limited streamer tubes for tracking (angular resolution 
~ 0.5°) and roughly 600 tons of liquid scintillator for fast timing (time resolution ~ 500 ps). The detector has overall 
dimensions of 12 x 77 x 9 m 3 and it is divided in 6 parts called supermodules. The bottom part of the apparatus, 
4.8 m high, is filled with rock absorber which sets a minimum threshold of about 1 GeV for vertical muons crossing 
the detector. The upper part, called the attico, is an open volume containing electronics as well as active detector 
elements. The streamer tubes form 14 horizontal and 12 vertical planes and the liquid scintillator counters form 3 
horizontal and 4 vertical planes on the outer surfaces of the detector. 

Neutrino events in the MACRO detector are seen in three different topologies: 

1. through-going upward muons. These events are produced by neutrinos interacting in the rock below MACRO 
and pass entirely through the detector. For atmospheric neutrinos the spectrum ranges from 1 to 10 4 GeV and 
the peak energy is about 100 GeV; 

2. internally produced upward-going events induced by neutrinos interacting in the lower part of the detector and 
producing a lepton which moves upward through the two upper scintillator layers. For atmospheric neutrinos 
the peak energy is about 4 GeV. To detect these events the attico must be in operation; 
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3. externally produced neutrino-induced upward-going stopping muons and internally-produced downward-going 
leptons. For atmospheric neutrinos the peak energy is about 4 GeV. 



The first and second topologies are recognized using the time-of-flight technique. Topological criteria are used for 
the third category. The time-of-flight technique is used to discriminate upward-going neutrino-induced muons from 
the background of downward-going atmospheric muons. Each scintillator records the time of a particle crossing by 
measuring the mean time at which signals are observed at the two ends of the scintillator. The transit time of a particle 
is found by taking the difference between the crossing times of two scintillators. In our convention, upward-going 
muons are considered to have a negative velocity, and downward-going muons a positive velocity. 

We present the results of WIMP searches using data gathered from March 1989 to March 1998. This data set 
encompasses five years of running with partial apparatus (March 1989 to April 1994) and 4 years of running with 
the full apparatus including the attico (April 1994 to March 1998) and it corresponds to livetimes including efficiency 
of 1.38 yr of running with the I s * super-module, 0.41 yr of running with the lower part of the detector and 3.1 yr of 
running of the full detector. The details of the upward-going muon analysis for the I s * supermodule and the lower 
detector can be found in Ref. p2| , for the full apparatus in Ref. |23] . 

Slightly different data sets and selection criteria are used for WIMP searches for the Earth and Sun to optimize the 
signal over background for each. For the Earth analysis only through-going upward muons which traverse at least 200 
g/cm 2 of absorber are used (category 1 above). The absorber requirement reduces the background due to soft pions 
produced at large angles by undetected downward-going muons to 1% Q. A total of 517 through-going upward 
muon events are used in the search for WIMP annihilation neutrinos from the Earth. 

Background rejection is not so critical for moving sources as it is for steady sources, so the analysis for the Sun 
includes through-going upward muons with no absorber requirement as well as semi-contained events generated in 
the lower half of the detector (category 2 above), for a total of 762 upward-going muon events. Apart from this, the 
analysis for the Sun is basically the same as the one for the Earth. 

We consider several contributions to the background of through-going upward muons: events with an incorrect 
timing measurement (such as muons in coincidence with radioactivity, other muons, or electromagnetic showers) and 
soft pions produced at large angles. We estimate 20 background events in the Earth sample of 517 through-going 
upward muons. 

The expected background of upward-going muons from atmospheric neutrinos is calculated with a full Monte Carlo 
calculation (described in Ref. ]22| , p3| ] ) . This calculation uses the Bartol flux |25|], the Morfin and Tung parton set Si 
p6[ for the deep inelastic vN cross-section and the muon energy loss in the rock from Ref. |27[ . We estimate a total 



uncertainty in the calculation of 17% 22 2^J. We have considered scenarios both with no neutrino oscillations as well 
as with the neutrino oscillation parameters of Ref. |p3| . The oscillations parameters considered are Am 2 = 2.5 x 10~ 3 
eV 2 and sin 2 29 = 1. 

In the no oscillation scenario the expected number of atmospheric neutrino events is 662 ± 113theor as compared 
to 462 ± 79theor in the oscillation scenario. In both cases there is a deficit of the measured events with respect to the 
expected number in the region around the vertical direction of MACRO, where the efficiency and acceptance of the 
apparatus are best known. Because of the discrepancy between the observed and expected numbers of upward-going 
muons in the expected signal region (0° to 30° from the vertical), we normalize the expected signal using data outside 
the expected signal region. This normalization is motivated by the fact that the absolute error of the expected flux is 
relatively high, whereas the shape of the flux is known to a few percent [p8[. This normalization factor is determined 
separately for each of the search cones considered, using the ratio of observed to expected events outside the search 
cone. 

We show both the expected atmospheric neutrino background with and without oscillations in Fig. ^(a) and (b). If 
neutrinos oscillate with these parameters the expected number of events would be reduced and the angular distribution 
of through-going upward muons would be distorted because neutrinos at the nadir oscillate more than at the horizon 
due to the longer pathlength (for a full discussion see Ref. Q ) . We note that the slight excess observed in the region 
—0.7 < cos0 < —0.6 is very unlikely to be produced by any plausible WIMP model, since this region is far from the 
expected signal region (9 < 30° - see Fig. |](a)). 

Muon flux limits are evaluated as: 

$ (90%c.l.) = N ^ 9 °/ OC - L ) (3.1) 

where N-p is the upper Poissonian limit (90% c.l.) given the number of measured events and expected background 
p9[ due to atmospheric neutrinos and £ is the exposure given by equation: 

£ = f e(t) x A(Q(t))dt (3.2) 
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where A(f2(t)) is the detector area in the direction of the expected signal (fi) at time t; e is the detector efficiency 
(discussed in Ref. [p2 23 1) which takes into account the possible variations of detector running configuration during 
data taking; and T sta rt and T cnc i are the start and end times of data taking. In the case of the Earth, the signal is 
expected always from the same direction; hence f2 is a constant with time and the detector acceptance A is calculated 
using the Monte Carlo calculation in Ref. J2^,^3|; the live-time and efficiency are calculated using downward-going 
events. 

Since the number of detected events is less than expected from the atmospheric neutrino flux, we set conservative 
flux limits assuming that the number of measured events in the signal region equals the number of expected events in 
that region J2j|. In Fig. |l|(b) the 90% c.l. muon flux limits for the Earth are plotted as a function of the nadir angle 
for 10 search half-cones around the vertical from 3° to 30°. For this data, the average value of the exposure is 2620 
m 2 yr (it varies slightly for different search cones because the area decreases by 37 m 2 when the search cone increases 
from 3° to 30° around the vertical). Note that the flux limits are independent on any hypothesis on a WIMP signal 
(or any other source). As shown in Fig. [l](b), the application of a — v T oscillation hypothesis to the atmospheric 
neutrino background with Am 2 — 0.0025 eV 2 and maximum mixing will result in lower flux limits. 

In the analysis for muons pointing in the direction of the Sun, the data themselves are used to generate the expected 
events in order to properly include the effects of semi-contained events. The arrival times from measured downward- 
going muons from the entire period of data taking are assigned randomly to the local trajectory coordinates of the 
measured upward-going muon events to evaluate their right ascension. This procedure allows us to take into account 
drifts of detection efficiency in time. Fig. ^(a) shows the angular distributions of measured and expected upward-going 
muons with respect to the direction of the Sun. The shape depends on the seasonal variation of the position of the 
Sun and on the livetime of the apparatus. The upward-going events detected during the night fall towards a cosine 
of 1, while the events collected during the day fall near -1. The exposure of the Sun is calculated with eq. 3.2 using 
the detector acceptance from the Monte Carlo described in Ref. |^] (for > 2 GeV, see below). The dependence 
of the acceptance of the detector on the direction of the Sun when it is below the horizon is exactly calculated using 
both the Monte Carlo and downward-going muon data. 

In Fig. |(b) the flux limits (90% c.l.) for an exposure of ~ 890 m 2 yr (it varies slightly for different search cones) for 
10 search cones around the Sun direction are shown. Again, it should be noted that any neutrino oscillation effects 
are automatically included in our method of deriving the number of expected events from real data. 

Table | shows the number of detected and expected upward-going muon events and corresponding flux limits for 
several cone sizes for the Earth and Sun. The muon threshold energy is determined by the amount of absorber 
an upward-going muon must traverse in MACRO. It is lower for the Earth (where tracks are oriented towards the 
vertical) than it is for the Sun (where tracks are more inclined). The upper limits in Tab. | are calculated assuming a 
minimum energy of 1.5 GeV and 2 GeV for the Earth and the Sun, respectively. Given the analysis requirements to 
select upward-going muons above these energies, we make a maximum error of 5% with respect to an exact calculation 
which takes into account the dependence on energy of the acceptance of the apparatus and of the neutrino fluxes from 
neutralinos with different masses. 



IV. FLUX LIMITS COMPARED TO SUSY PREDICTIONS 



The flux limits from the previous section can be used to constrain any WIMP model. In this section, we apply our 
flux limits to the neutralino models of Bottino et al. |l(J using search cones which collect 90% of the expected signal. 
The range of search cones accommodates a range of neutralino masses. 

The angular distribution of upward-going muons follows that of the parent neutrinos, with deviations due to charged- 
current neutrino interaction and Coulomb multiple scattering of the produced muons in their path to the detector. 
Both of these distributions are a function of the neutrino energy. In turn, the energy spectrum of neutrinos produced 
in WIMP annihilation depends both on the WIMP mass and the annihilation products. The final states can be 
pairs of fermions, or Higgs/gauge bosons or combinations of Higgs and gauge bosons. The branching ratios into these 
channels depend on the model and they have some effect on the neutrino energy spectra. Annihilation to fermion pairs 
tends to produce softer neutrinos than annihilation to gauge bosons since fermions dissipate energy in hadronization 
whereas bosons have a greater likelihood of prompt decay to neutrinos. We estimate the maximum variation in flux 
limits by taking the extreme cases where only one annihilation channel is open. We find this variation is no more 
than 17% between extreme models and those considered here. To a good approximation, the angular distribution is 
a function of just the neutralino mass, and so we determine flux limits as a function of neutralino mass. 

The angular distributions of the upward-going muon signals are calculated using neutrino fluxes from neutralino 
annihilation in the Sun and Earth calculated by Bottino et al. [flQf| . A Monte Carlo calculation with cross-sections 
described in Ref. 1 30 1 and muon energy loss as described in Ref. plj] is used to propagate muons through the detector 
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where the angular resolution is taken into account. 

The main difference between the angular distributions of the signals from Earth and Sun is due to the angular 
resolution of the detector which degrades for slanted tracks because the number of streamer tube layers crossed 
decreases. Hence, the more slanted tracks from the Sun tend to have lower angular resolution than those from the 
Earth. Moreover, the angular size of the expected signal is affected by the angular size of the region where WIMP 
annihilation is taking place. The diameter of the Sun is 0.5° (as seen by the Earth) and the WIMPs are expected 
to be localized at the center of the Sun. As a result, the angular size of the WIMP annihilation region is negligib le. 
However, for the Earth the angular size of the annihilation region is considerable and has been estimated by J32|,|l7],[l5| 
to be: 

G(6) ~ 4m x ae~ 2m * Qsin2 9 (4.1) 

where 9 is the nadir angle; a is a parameter depending on the central temperature (T = 6000 K), the central density 
(p = 13 g cm" 3 ) and radius of the Earth (a = 1.76 GeV" 1 ). 

In Fig. |||the muon nadir angle is shown for m x — 60, 100, 200, 500 and 1000 GeV and fixed model parameters. In 
Fig. |^ the angular spreads between the neutrino and the muon directions are shown as a function of neutralino mass 
in the case of the search for neutralinos trapped inside the Sun. 

The 90% c.l. flux limits are calculated as a function of neutralino mass using cones which collect 90% of the 
expected signal. These limits are corrected for the 90% collection efficiency due to cone size. Figures || and ^ show 
these limits for the Earth and Sun, respectively. The experimental limits are superimposed on the flux of upward-going 
muons from the Bottino et al. calculation as a function of m x . The flux limits are calculated for the same muon 
minimum energy as is used in the calculation (E^ > 1 GeV) considering the dependence on the energy of the MACRO 
acceptance in the low energy region. A correction is applied for each neutralino mass to translate from the thresholds 
of Tab. | to the 1 GeV threshold used in the calculation of upward-going muon fluxes. The correction factors are 
higher for lower neutralino masses. Moreover, they are higher for the Sun for which the threshold in Tab. | is 2 GeV 
than for the Earth for which it is 1.5 GeV. We estimate these factors to be 1% for the Earth and 10% for the Sun for 
a neutralino mass of 60 GeV. 

The fluxes are calculated by varying the model parameters (each model is represented by a dot in Fig. || and in 
Fig. ^ in experimentally allowed ranges. Dots correspond to variations of parameters between the following values: 
tan/3 = 1.01, 2, 3, 10, 40, 50; m A = 65, 70, 75, 80, 85, 90, 95, 100, 200, 300, 500 GeV, where m A is the mass of the 
pseudoscalar Higgs; m = 150, 200, 300, 500 GeV, where mo is the common soft mass of all the sfermions and squarks; 
A = -3, -1.5, 0, 1.5, 3 where A is the common value of the trilinear coupling in the superpotential for the bottom and 
top quark (A is set to zero for the first and second family); \p,\ and Mi are varied between 10 and 500 GeV in steps of 
20 GeV. The calculation also assumes a neutralino rms halo velocity of 270 km s _1 , a halo escape velocity of 650 km 
s -1 , a velocity of the Sun around the Galactic center of 232 km s _1 , a local dark matter density of 0.5 GeV cm" 3 , 
and a minimal value for rescaling the neutralino relic abundance, (ilh 2 ) m i n , of 0.03 (explained below). Those fluxes 
lying above the experimental flux limit curve are ruled out as possible SUSY models by this measurement given the 
cosmological parameters chosen. It should be considered that a variation of the astrophysical parameters may lower 
the calculated fluxes by at most one order of mag nitude Moreover, — v r oscillations could lower neutrino 

fluxes from neutralino annihilations by about a factor of two for the oscillation parameters considered in this paper 

& 1-1 

MACRO flux limits can be used to constrain also other neutralino calculations, e.g. Bergstrom et al. |16fl. In this 

calculation configurations with Q x h 2 < 0.025 are not taken into account. On the other hand, Bottino et al. ||l0||, when 
the neutralino cosmological abundance is too low to account for the total dark matter in the halo, assume: 



if tt x h 2 > {Wi 2 ) min p x = pioc 

if VL x h 2 < (Wi 2 ) min p x = pioc x n x h 2 /(nh 2 ) r 



(4.2) 



where (Vlh 2 ) m i n — 0.03 and pi oc — 0.5 GeV cm" 3 . Both calculations indicate that indirect searches can have good 
prospects particularly at high neutralino masses (see Fig. |^ and Fig. ||). 

We have tested another approach to increase the sensitivity in the estimate of flux limits for the Sun through the 
study of the angular distribution of the signal |34| . The data, signal, and background angular distributions are fit in 
the region of interest for the signal, using a x 2 expression suited for Poisson-distributed data j29| . In this expression 
the signal is multiplied by a factor K . The normalization factor K of the signal and its error are found by minimizing 
X 2 with respect to K . The resulting flux limits are evaluated at the 90% c.l. as a function of the neutralino mass. 
The results of this method, although still with limited statistics, are in agreement with the previous one in which the 
flux limits are calculated in cones which collect 90% of the signal, because our data in the direction of the Sun are in 
good agreement with the prediction of the atmospheric neutrino background. In principle, this method is eventually 
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preferable because it avoids an a priori choice of the search cone and hence of the fraction of the signal lost outside 
of it and it will improve with increasing statistics. 

Recently the DAMA/Nal experiment observes a possible a possible annual modulation effect in a WIMP direct 
search at a 99.6% c.l. |35[| . This modulation has been interpreted in terms of a relic neutralino which may make up the 
major part of dark matter in the universe (see Ref. |3(| and references therein). Fig. [7] shows the allowed SUSY models 
considered by Bottino et al. for various neutralino local density compared to the MACRO upward-going muon flux 
limits from the Earth pq] . MACRO experimental upper limits from the Earth rule out many SUSY configurations 
indicated by the DAMA/Nal experiment, even assuming the atmospheric neutrino background oscillates with the 
parameters favored by MACRO [^3| and Superkamiokande p7[ . The MACRO data from the Sun, which is more 
sensitive to spin-dependent scattering, have less overlap in sensitivity with direct searches. We again notice that 
expected fluxes may decrease by about a factor of two in the presence of neutrino oscillations. 



V. CONCLUSIONS 



A search for a WIMP signal is performed using the MACRO detector at the Gran Sasso Laboratory to observe 
upward-going muons coming from the products of WIMP annihilations from the Earth and Sun. The time-of-flight 
technique is used to discriminate the upward-going neutrino-induced events from the background of downward-going 
atmospheric muons. 

We look for an excess of neutrino events over the background of atmospheric neutrinos in the direction of the Sun 
and Earth. No signal of WIMP annihilation is observed, and we present upper limits on the flux of upward-going 
muons from these bodies. 

While we do not see a WIMP signal, we compare our flux limits to predictions from various SUSY models and rule 
out those which are inconsistent with our limits. We calculate the expected upward-going muon fluxes and angular 
distributions using the SUSY neutralino annihilation calculations of Bottino et al. p0| . The 90% c.l. upper limits for 
the flux of upward-going muons in angular cones collecting 90% of the expected signal from neutralino annihilation 
are given. Our data exclude significant portions of the parameter space for neutralinos and our limits for the Earth are 
the most stringent of all "indirect" experiments 1 18 ljj . The flux limits for the Sun however tend to rule out a smaller 
portion of the SUSY parameter space at this time. Finally, MACRO data also exclude some of the SUSY models 
suggested by the annual modulation analysis performed by the DAMA/Nal experiment (especially at lower neutralino 
masses) p5fl . The MACRO indirect search for WIMPs has reached the proper sensitivity to explore configurations 
compatible with the DAMA/Nal data. 



ACKNOWLEDGMENTS 



We thank Prof. A. Bottino and Dott. N. Fornengo for valuable discussions and for their contributions to this 
work. We gratefully acknowledge the staff of the Laboratori Nazionali del Gran Sasso and the invaluable assistance 
of the technical staffs of all the participating Institutions. For generous financial contributions we thank the U.S. 
Department of Energy, the National Science Foundation, and the Italian Istituto Nazionale di Fisica Nucleare, both 
for direct support and for FAI grants awarded to non-italian MACRO collaborators. 



Corresponding author E-mail: montaruli@ba.infn.it 
[1] J. R. Primack, in "Dark Matter in Astro- and Particle Physics", edited by H. V. Klapdor-Kleingrothaus and Y. Ramachers 

(World Scientific, Singapore, 1997), p. 97. 
[2] A. Dekel and M. J. Rees, Astrophys. J. 422, LI (1994). 
[3] K. A. Olive, in "Dark Matter in Astro- and Particle Physics" Q, p. 250. 

[4] G. F. Smoot et al., Astrophys. J. 396, LI (1992); C. L. Bennett et al, Astrophys. J. 464, LI (1996). 
[5] M. Persic, P. Salucci, and F. Stel, Mon. Not. Roy. Astron. Soc. 281, 27 (1996). 
[6] J. P. Ostriker and P. J. E. Peebles, Astrophys. J. 186, 467 (1973). 
[7] C. Alcock et al, Astrophys. J. 486, 697 (1997). 

[8] G. Jungman, M. Kamionkowski, and K. Griest, Phys. Rep. 267, 195 (1996). 
[9] H. E. Haber and G. L. Kane, Phys. Rep. 117, 75 (1985). 



7 



A. Bottino, N. Fornengo, F. Donato, and S. Scopel (private communication). N. Fornengo, to appear in the Proceedings 
of the Ringberg Euroconference "New Trends in Neutrino Physics" (Ringberg Castle, Tegernsee, Germany, May 1998), 
edited by B. Kniel (World Scientific, Singapore). 

OPAL Collaboration, K. Ackerstaff et al, Eur. Phys. J C 2, 213 (1998); ALEPH Collaboration, R. Barate et al, Eur. 

Phys. J. C 2, 417 (1998). 

K. Freese, Phys. Lett. B 167, 295 (1986). 

L. Krauss, M. Srednicki, and F. Wilczek, Phys. Rev. D 33, 2079 (1986). 

J. Silk, K. A. Olive, and M. Srednicki, Phys. Rev. Lett. 55, 257 (1985). 

A. Bottino, N. Fornengo, G. Mignola, and L. Moscoso, Astropart. Phys. 3, 65 (1995). 

L. Bergstrom, J. Edsjo, and P. Gondolo, Phys. Rev. D 55, 1765 (1997). 

A. Gould, Astrophys. J. 321, 571 (1987). 

BAKSAN Collaboration, M. M. Boliev et al, Nucl. Phys. B 48 (Proc. Suppl.), 83 (1996); and in "Dark Matter in Astro- 
and Particle Physics" M, p. 706. 

KAMIOKANDE Collaboration, M. Mori et al, Phys. Rev. D 48, 5505 (1993). 
1MB Collaboration, J. M. LoSecco et al, Phys. Lett. B 188, 388 (1987). 
MACRO Collaboration, S. P. Ahlen et al, Nucl. Inst, and Meth. A 324, 337 (1993). 
MACRO Collaboration, S. P. Ahlen et al, Phys. Lett. B 357, 481 (1995). 
MACRO Collaboration, M. Ambrosio et al, Phys. Lett. B 434, 451 (1998). 
MACRO Collaboration, M. Ambrosio et al, Astropart. Phys. 9, 105 (1998). 
V. Agrawal, T. K. Gaisser, P. Lipari, and T. Stanev, Phys. Rev. D 53, 1314 (1996). 
J. G. Morfin and W. K. Tung, Z. Phys. C 52, 13 (1991). 

W. Lohmann et al, CERN Yellow Report No. CERN-EP/85-03, 1985 (unpublished). 

P. Lipari and M. Lusignoli, Phys. Rev. D 57, 3842 (1998). 

C. Caso et al, Review of Particle Physics, Eur. Phys. J. C 3, 1 (1998). 

P. Lipari, M. Lusignoli, and F. Sartogo, Phys. Rev. Lett. 74, 4384 (1995). 

P. Lipari and T. Stanev, Phys. Rev. D 44, 3543 (1991). 

K. Griest and D. Seckel, Nucl. Phys. B 283, 681 (1987). 

N. Fornengo (private communication). 

J. Edsjo and P. Gondolo, Phys. Lett B 357 (1995) 595; L. Bergstrom et al, Astropart. Phys. 7, 147 (1997). 
R. Bernabei et al, Phys. Lett. B 424, 195 (1998); and University of Rome report No. ROM2F/98/34 and INFN/AE-98/20 
(1998). 



A. Bottino, N. Fornengo, F. Donato, and S. Scopel, hep-ph/9809239, to appear in Astropart. Phys.. 
SUPERKAMIOKANDE Collaboration, Y. Fukuda et al, Phys. Rev. Lett. 81, 1562 (1998). 
R. Bernabei et al, Phys. Lett. B 389, 757 (1996). 

BESS Collaboration, H. Matsunaga et al, in the Proceedings of the 25 th International Conference of Cosmic Rays (30 July 
- 6 August 1997, Durban, South Africa), edited by M. S Potgieter, B. C. Raubenheimer and D. J. van der Walt, (Space 
Research Unit, 1997) vol. 3, OG Sessions 1-5, p. 373. 



EARTH SUN 



Cone 


Data 


Back- 
ground 
events 


Norm, 
factor 


Flux Limit 
(£ M > 1.5 GeV) 
(cm" 2 s" 1 ) 


Data 


Back- 
ground 
events 


Flux Limit 

(E„ > 2 GeV) 
(cm" 2 s" 1 ) 


30° 


76 


119.2 


0.82 


2.24 xl0~ 14 


56 


51.8 


5.84 xlO" 14 


24° 


52 


75.5 


0.79 


1.75 xlO" 14 


33 


33.1 


3.83 xlO" 14 


18° 


32 


42.2 


0.77 


1.42 xl0~ 14 


17 


18.8 


2.60 xlO -14 


15° 


24 


29.0 


0.76 


1.20 xlO" 14 


11 


13.2 


2.09 xl0~ 14 


9° 


10 


10.5 


0.75 


7.48 xl0~ 15 


3 


4.7 


1.34 xl0~ 14 


6° 


4 


4.6 


0.75 


5.30 xlO -15 


2 


2.1 


1.37 xlO" 14 


3° 





1.1 


0.75 


3.79 xlO -15 


2 


0.5 


1.73 xl0~ 14 



TABLE I. Observed and background events and 90% c.l. muon flux limits for some of the 10 half-cones chosen pointing 
from the Earth and the Sun. The background events are those expected from atmospheric neutrinos. In the case of the Earth, 
the normalization factors are used to normalize the expected background events by the ratio of observed to expected events 
outside each cone. Since the number of detected events (column 2) is less than the normalized expected events (column 3), 
we set conservative flux limits assuming that the number of measured events equals the number of expected ones The 
Earth results are for the no oscillation scenario. Earth limits considering neutrino oscillations agree to within 8%. The average 
exposure for the Earth is 2620 m 2 yr and for the Sun 890 m 2 yr. 
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FIG. 1. (a) Nadir distribution of measured (black circles) and expected (solid line and dashed line) through-going upward 
muons. The expected distributions are multiplied by the ratio of the measured events over the expected ones outside the 
largest window (30°). The normalization factors are 0.82 for the solid line (no oscillations) and 1.19 for the dashed line (y^ — v T 
oscillations for maximum mixing and Am 2 = 0.0025 eV 2 ). (b) Muon flux limits (90% c.l.) as a function of the angle from the 
vertical (the angle varies form 3° to 30° in steps of 3°). In both plots, the dashed line is obtained in the hypothesis of — > v T 
oscillations of the atmospheric neutrino background with maximum mixing and Am 2 = 0.0025 eV 2 . 
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FIG. 2. (a) Distribution of measured (black circles) and expected (solid line) upward-going muons as a function of the cosine 
of the angle from the Sun direction, (b) Muon flux limits (90% c.l.) as a function of the search cone around the direction of 
the Sun. 
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FIG. 3. Nadir angle distribution of muons induced by neutralino annihilation inside the Earth for several neutralino masses. 
The angular ranges including 90% of the signal are indicated. 
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FIG. 4. Neutrino-muon angular separation distribution for neutrinos from x~X annihilation in the Sun for several neutralino 
masses. The angular ranges with 90% of the signal are shown. 
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FIG. 5. Upward-going muon flux vs m x for _E* h = 1 GeV from the Earth Jl^. Each dot is obtained varying model parameters. 
In this plot values of fi > are considered. Similar results are obtained for fi < 0. Solid line: MACRO flux limit (90% c.L). 
The solid line representing the flux limit for the no-oscillation hypothesis is indistinguishable in the log scale from the one for 
the Vfj, — v T oscillation hypothesis, but the expectations could be about two times lower. The open circles indicate the models 
excluded by direct measurements (particularly the DAMA/Nal experiment S|) and assume a local dark matter density of 
0.5 GeV cm -3 . See Fig. |^ for the comparison for the same density between the MACRO flux limit and the allowed values of 
parameters based on recent DAMA/Nal results. 



12 




m x (GeV) 
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Solid line: MACRO flux limit (90% c.L). The open circles concern the regions excluded by direct measurements 
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FIG. 7. The dashed lines are our experimental upgoing muon flux limits (for E 1 ^ = 1 GeV) at 90% c.l. as a function of 
m x for the Earth. The dots and the circles are SUSY neutralino models allowed at 90% c.l. from Bottino et al. assuming 
various local dark matter densities J3^] . The experimental flux limit and the theoretical fluxes shown are obtained assuming no 
neutrino oscillations. The limits computed assuming oscillations are indistinguishable in this graphs; the expected fluxes may 
decrease by about a factor of two. The dots represent models already excluded by cosmic anti-proton data [M . 
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